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Background: The process of /3-delayed fission (/3DF) provides a versatile tool to study low-energy fission in 
nuclei far away from the /1-stability line, especially for nuclei which do not fission spontaneously. 

Purpose: The aim of this paper is to investigate systematic trends in /IDF partial half-lives. 

Method: A semi-phenomenological framework was developed to systematically account for the behavior of /IDF 
partial half-lives. 

Results: The /IDF partial half-life appears to exponentially depend on the difference between the Q value for 
j3 decay of the parent nucleus and the fission-barrier energy of the daughter (after /3 decay) product. Such 
dependence was found to arise naturally from some simple theoretical considerations. 

Conclusions: This systematic trend was confirmed for experimental /3DF partial half-lives spanning over 7 orders 
of magnitudes when using fission barriers calculated from either the Thomas-Fermi or the liquid-drop fission model. 
The same dependence was also observed, although less pronounced, when comparing to fission barriers from the 
finite-range liquid-drop model or the Thomas-Fermi plus Strutinsky Integral method. 
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Figure 1. (Color online) Schematic representation of the /IDF process on the neutron-deficient side of the nuclear chart. 
The Qec value of the parent (A,Z) nucleus is indicated, while the curved line shows the potential energy of the daughter 
(A,Z-1) nucleus with respect to nuclear elongation, displaying also the fission barrier Bf. The color code on the right-hand 
side represents the probability for excited states, with excitation energies close to Bf, to undergo fission; the darker colors 
correspond to higher probabilities. 


I. INTRODUCTION 

,0-delayed fission (/3DF) is a two-step process whereby the fissioning nucleus could be created in an excited state 
after f3 decay of a precursor. Since the excitation energy of the fissioning daughter product is limited by the Qjs 
value for /? decay of the parent, /3DF provides a unique tool to study low-energy fission of nuclei far from stability, 
especially for those not fissioning spontaneously. Figure provides a schematic representation of this process, for 
nuclides on the neutron-deficient side of the nuclear chart. Recent experiments at ISOLDE-CERN [IHl] and SHIP- 
GSI [3 [6] have studied this exotic decay mode in several short-lived neutron-deficient isotopes in the lead region. 
The fission-fragment mass and energy distributions resulting from /3DF have established a new region of asymmetric 
fission around mm and indicated multimodal fission in i94,i96pQ ^^^Rn |3]. A recent review of the /3DE 

process is given in [7], in which a total of 27 /3DF cases, both on the neutron-rich and neutron-deficient sides, were 
summarized. 

It is furthermore believed that /3DF could, together with neutron-induced and spontaneous fission, influence the 
fission-recycling in r-process nucleosynthesis Therefore, a reliable prediction of the relative importance of ,5DE 

in nuclear decay, often expressed by the /3DF probability P/ 3 df, is needed. P/jdf is defined as 


P/3DF 


NjSDF 


( 1 ) 


where Npup and A(g are respectively the number of /3DF and j5 decays of the precursor nucleus. An earlier 
comparison of P/jdf data in a relatively narrow region of nuclei in the vicinity of uranium showed a simple exponential 
dependence with respect io Qp [TOllII]- R was assumed that fission-barrier heights Bf of the daughter nuclei do not 
vary greatly in this region [T2] {Bf ~ 4-6MeV) and thus have a smaller influence on PpuF as compared to Qp values 
{Qp ~ 3-6MeV). In addition, these nuclei have a typical N/Z ratio around ^1.4-1.5, which is close to that of 
traditional spontaneous fission of heavy actinides. 

The aim of this paper is to further explore such systematic features by including the newly obtained data in the neutron- 
deficient lead region whose /?DF nuclides have significantly different N/Z ratios (^1.2-1.3), Bf (~7-10MeV) and 
Qp values (^9-IlMeV) as compared to those in the uranium region. 

However, from an experimental point of view, the dominant a-branching ratio (> 90 %) in most /3DE precursors in the 
neutron-deficient lead region |13j makes precise determination of Np in equation Q difficult. Therefore, the partial 
/3DF half-life T’i/ 2 p,/ 3 DF) as proposed in [7], is discussed in the present study. By analogy with other decay modes, 
Ti/ 2 p,pF>F is defined by 


7’l/2p,;SDF 


NpDF 


( 2 ) 


where T 1/2 represents the total half-life and A/jec.tot the number of decayed precursor nuclei. The relation between 
7’i/2p,/3DF and P^df can be derived from equations Q and ([^ as 
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^l/2p,/3DF 


T"i/2 

bjjPp-DY ’ 


(3) 


with bp denoting the /3-branching ratio. If the a-decay channel dominates, as is often the case in the neutron- 
deficient lead region, one can safely approximate IVdec.tot in equation ® by the amount of a decays 
This work shows an apparent exponential dependence of 7 i/ 2 p,/ 3 DF on — B{) for certain sets of calculated fission- 
barrier energies. Such relation may arise naturally by simple phenomenological approximations of the /3-strength 
function of the precursor and the fission-decay width of excited states in the daughter nucleus. These assumptions 
may be justified considering the scale of the systematic trend discussed here, spanning Ii/ 2 p„SDF values over several 
orders of magnitude. Deviations lower than one order of magnitude are thus acceptable. 


II. THEORETICAL CONSIDERATIONS 
Following [HHis], the expression for PpuF is given by 


_ CSpiE)F{Qp-E)^dE 

Sp{E)F{Qp-E)dE 

whereby the /3-strength function of the parent nucleus is denoted by Sp and the Fermi function by F. The 
excitation energy is here, and further, given by E. The fission and total decay widths of the daughter, after /3 decay, 
are respectively given by Ff and Ftot- Equation ([^ can be combined with equation Q to deduce the decay constant 
of /3DF, defined as A/jdf = ln( 2 )/Ti/ 2 p,/ 3 DF, as 


Tt(E) 

XpDF= SpiE)F{Qp-E)-^^dE. (5) 

Jo Itot(Ej 

This section will be devoted to the derivation of an analytical expression for A^df by approximating Sp, F and 
Ff/Ftot- Since most of the reliable experimental data on /3DF are recorded on the neutron-deficient side of the nuclear 
chart (see Table and H), only EC//3“''-delayed fission will be considered here. 


A. Approximations 

A first simplification in equation ([^ is to approximate Sp by a constant Ci, as proposed in previous studies (see 
for example nzmsi). Possible resonant structures in Sp, considered in e.g. [151 IIS]j are thus ignored, thereby 
assuming a limited sensitivity of Ti/ 2 p,/ 3 DF on Sp with respect to the scale of the systematic trend discussed here. 
This approximation is further supported by the study in |20j . which shows a limited influence of Sp in the calculation 
of PpDF- Furthermore, Ci was taken equal for all isotopes listed in Table [Tj thereby neglecting possible variations of 
Cl with respect to the nuclear properties of the /3DF precursors - such as mass, proton number, isospin, spin and 
parity. 

The Fermi function F can be fairly well described by the function C 2 {Qec ~ PY [SHIS] for EC decay. The prefactor 
C 2 was again considered element independent, thereby ignoring its slight dependence on the atomic number Z [23j . 
According to [551 [23], EC decay is dominant for transition energies below 5MeV if Z exceeds 80. Since Qp values of 
/3DF precursors in the uranium region are typically smaller than 5MeV (see Table [^, decay can be disregarded 
here. Qp values in the neutron-deficient lead region can however reach 10-llMeV, implying a relatively high 
over EC decay ratio to the ground or a low-lying excited state in the daughter. However, since /3DF should primarily 
happen at excitation energies which are only a few MeV below Qp [55], EC-delayed fission should dominate over /3+ 
delayed fission in the full region of the nuclear chart (see further). 

The prompt decay of an excited state in a nucleus can, in the most general case, happen through fission, emission of 
a 7 ray, proton, a particle or neutron. The total decay width is thus given by Ftot = Ff -|- F..^ -|- Fp -|- Fa -|- Fn. 

For the /3DF precursors considered in Table [Tj the neutron separation energies exceed the Qp value by at least several 
MeV [26] and charge particle-emission is strongly hindered due to the large Coulomb barrier. Therefore, the de¬ 
excitation of states below is mostly dominated by 7 decay, which makes that Ftot — r..y [5D1[57]. In addition, F.^ 
can be approximated by a constant (see for example [50]). Reference [57] provides a calculation of F/ with respect to 
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Figure 2. Plot showing the integrand of equation ([^, which is proportional to the /3DF probability, for X equal to 3,4 or 5. 


the excitation energy E by including the fission-barrier penetrability and the influence of level densities at the ground 
state and saddle point. This calculation shows that Tf seems well approximated by a single exponential behavior 
Ff ~ g-x{B!-E) excitation energies around B{. For the fissile nuclei listed in Table|^ the decay constant adopts a 
value X Ki A MeV“^ [57]. The ratio Ff/Ttot is thus approximated by 


^ ( 6 ) 

-L tot ^ 7 

The constants and X are assumed to adopt the same value for all isotopes of interest. At excitation energies E 
moderately above Bf, de-excitation by fission should dominate and Ti/Ttot{E) will thus be close to unity. Since the 
Qp value of most known /3DF precursors (see Table 1^ does not exceed B{ of the daughter by more than a few MeV, 
it is further assumed that equation remains valia for excitation energies in the daughter nucleus close to 
Using the above approximations ancftaking C = C 1 C 2 C 3 , the right-hand side of equation ([^ reduces to 

rQfi 

A;3df = C (Qp- (7) 

Jo 


B. Calculating Xtdf 

Equation Q can be rewritten, in order to isolate the exponential dependance on (Qp — i?f), as 

rQfi 

XpYiF = / [Qp - Efe-^^^f^-^'^dE. ( 8 ) 

^0 

The integrand in equation Q is thus proportional to the /3DF probability at a given E of the daughter nucleus. This 
function, plotted in Figure]^ for different values of X around the deduced value X « 4MeV“^ from |27|, shows that 
/3DF primarily happens at energy levels 0-2MeV below Qp. Moreover, since all Qp values of the neutron-deficient 
/3DF precursors listed in Table |I] exceed ~ 2 MeV, the value of the integral in equation Q is little dependent on the 
precise value of Qp. As a consequence, A^df primarily depends on the difference [Qp — B[). 

In order to prove latter statement analytically, a substitution with u = X[Qp — E) and adjustment of integration 
borders in equation ([^ is performed: 


A/3DF = 


rXQp 

./n 


u^e "^du. 


(9) 


The integral in equation (§ is similar to the mathematical form of the so-called normalized upper incomplete 
Gamma function, defined as 
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Figure 3. The normalized incomplete Gamma function V(Z, XQp), needed for the calculation of the integral under the /3DF 
probability curves shown in Figure 


r(s,x) 




r(s) 


dt, 


whereby r(s) is 


Equation Q thus transforms into 


f+oo 

r(s) = / ; 

Jo 


*dt. 


( 10 ) 


( 11 ) 


(JpX{Qi3-B!) 

X/3DF = — ^3 -r(3)r(3,xQ^). 


( 12 ) 


Table |T] shows that all Qp values of the neutron-deficient /3DF precursors exceed 3MeV, while the fitted values for 
X in Table [ill as well as the theoretical estimate from [27] {X « 4MeV“^), are all greater than 1.7MeV“^. The 
value XQp thus exceeds 5 in all discussed cases, implying that, as shown in Figure]^ one can thus safely approximate 
r(3,X(5^) ~ 1 in equation (12). 

In this simple picture, it isThus found that h3{XpDF) depends linearly on (Q^ — Bi). In terms of the partial /3DF 
half-life 7i/2p,/3DF one finds the relation 


loglo(T’i/2p,/3DF) — C — Xlogio(e)((5/3 — Bi), 


(13) 


with the constant C given by 


III. SYSTEMATIC COMPARISON OF EXPERIMENTAL DATA 

This section aims at verifying equation © by using experimental /3DF partial half-lives and theoretical values for 
(Q^ — B{), summarized in Table and Figure]^ Tabulated fission barriers from four different fission models were 
used, of which three are based on a macroscopic-microscopic and one a mean-field approach. The latter model is 
based on the Extended Thomas-Fermi plus Strutinsky Integral (ETFSI) method [55|, but tabulated barriers for the 
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Table I. List of all precursors for which /3DF was observed. The measured half-life T 1 / 2 , /3-branching ratio hp, /3DF probability 
PpuF, ratio of observed /3DF to a decays NpoF/Na and calculated /3DF partial half lives Ti/ 2 p,/ 3 df are listed. Reliable values 
for ri/ 2 p,/ 3 DF, as evaluated by the criteria in 0, are indicated in bold. (Qp — Bf) is tabulated for hssion barriers from four 
different fission models : Thomas-Fermi (TF) [^, Finite Range Liquid Drop (FRLDM) [ST], Liquid Drop (LDM) [27] and 
the Extended Thomas-Fermi plus Strutinsky Integral (ETFSI) model |28| . Qp values were taken from |5B] and are defined by 
equation (15l. 






Qp - Bt (MeV) 
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^ Value extracted according to equation ( |16[ l by using evaluated experimental data from m- 
Calculated /3-branching ratio from m- _ 

^ Value extracted according to equation ( |16[ | by using experimental data from ESj. 

^ /3-branching ratio from m- 
® Evaluated ^-branching ratio from m- 


most neutron-deficient isotopes in Table |T] are absent in literature. The microscopic-macroscopic approaches all rely 
on shell corrections from [29] and describe the macroscopic structure of the nucleus by either a Thomas-Fermi (TF) 
m. liquid-drop (LDM) or the Finite-Range Liquid-Drop Model (FRLDM) |3T]. The Qp values were taken from 
the 2012 atomic mass evaluation tables [26] and are derived from the difference between the atomic masses of parent 
Mp{Z,A) and daughter Mo{Z',A) nuclei as 


Qp = ^[Mp{Z,A)-Md{Z',A)\. (15) 

About half of these values are known from experiments, while the others are deduced from extrapolated atomic 
masses. In latter cases, the difference of the Qp values from [26] with the theoretical values from [29] or [32] is always 
lower than 0.4 MeV. 

2 \/ 2 p„flDF values were extracted from reported Ppp,p values using equation (|^, if the precursor nucleus has a 
significant /3-decay branch (bp > 10%). When multiple measurements on Ppjpp were performed, only the reliable 
value, as evaluated by [7], or the most recent value was tabulated. In case of a dominant a-decay branch (bp < 10 %), 
Ti/ 2 p,pDF was calculated by equation ([^, whereby Vdec,tot was approximated by the observed amount of a decays 
Nd, corrected for detection efficiency. 

Since the isotopes 192 , 194 ^^ 202 pj. have both a ground and a low-lying alpha-decaying isomeric state 
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Table II. Results of the fits, corresponding to four different fission models, shown in Figure]^ The values for the parameters X 
and C' in equation (13l are listed. Also the root-mean-square deviations (RMSD) of the reliable experimental logio(7i/2p,/3DF) 
values (represented by the closed symbols in Figure]^ to the fit are given. 


Model 

X (MeV-i) 

C’ (MeV) 

RMSD 

TF 

3.0(2) 

6.2(1) 

0.47 

FRLDM 

1.7(4) 

4.9(3) 

1.19 

ETFSI 

2.1(7) 

5.0(6) 

1.10 

LDM 

2.2(2) 

5.8(2) 

0.62 


with comparable half-lives, only an overall Nf^D-p/Na value could be extracted with present experimental techniques. 
We refer the reader for a detailed discussion of this issue to HE]. Therefore, these precursors have been excluded 
from the fit in Figure |4| Nonetheless, as a first approximation the value for Tx/ 2 p,/ 3 DF extracted by defining the 
half-life Ti/ 2 ^ shown in~table|^ as the unweighted average 


Ti/2 — 


Tl/2,g, + 7l/2,m 
2 


(16) 


where the respective half-lives for ground and isomeric states are denoted by T\i 2 ,^ and rL/ 2 ,m- The uncertainty 
Ari /2 is conservatively taken as 


^Ti/2 — 


|T’l/2,g — ^1/2,ml 
■ ^2 


(17) 


Figure 1^ shows logio('ri/ 2 p,/ 3 DF) against {Q^ — B{) for the fission barriers from the four different models under 
consideration. Using the same evaluation criteria as proposed in [7] for P/jdf measurements, 13 reliable '7i/2p,/3DF 
values, marked in bold in Table |l| were selected. These data points, represented by the closed symbols, are fitted by a 
linear function. An equal weight to all fit points is given because the experimental uncertainties on logic (7i/2p,/3df) 
are in most cases much smaller than the deviation of the data points with the fitted line, of which the extracted 
parameters are summarized in Table |TT| The remaining data points from Table |lj are shown by open symbols and were 
excluded from the fit. The color code discriminates between the neutron-deficient lead region (red), neutron-deficient 
(black) and neutron-rich (blue) uranium region. 

Figure lil illustrates a linear dependence of logic(T i/ 2 p,/ 3 df) on (Qp — Bf) for TF and LDM barriers for over 7 orders 
of magnitude of T'i/ 2 p,/ 3 DF- In addition. Table [ll| shows a relatively small root-mean-square deviation (RMSD) of the 
13 reliable experimental logic(T i/ 2 p,^df) values (represented by the closed symbols in Figure]^ to the corresponding 
values extracted from the fit. The dependence is somewhat less pronounced for the FRLDM model, as evidenced by a 
larger RMSD value. A similar linear trend is observed for the ETFSI model, but the lack of tabulated fission barriers 
in the neutron-deficient region, especially in the lead region, prohibits drawing definite conclusions. 

Moreover, Table [ll| shows that the four fitted values for X are similar to each other as well as to the theoretical 
estimate X « 4MeV“^ m- The extracted values for the offset C are also found to be comparable. 


In contrast to a rather good agreement for most neutron-deficient nuclei, all models show a larger systematical 
deviation from this linear trend for the neutron-rich /3DF precursors ^^®Ac and In [7], concerns were 

raised on the accuracy of the P^sdf values measured in this region, which could explain this deviation. Note also 
that the precursors in this region of the nuclear chart undergo j3~ decay in contrast to the EC-delayed fission on 
the neutron-deficient side for which equation (131 was deduced, influencing the numeric value of the offset C. In 
particular, the Fermi function for /3“ decay is approximately proportional to {Qp — E)^ [THIIHI, in contrast to the 
quadratic dependence on (Q^ — E) for EC decay. The parameter X should however remain unchanged, because 
equation (|^ approximating Ff/Ftot remains valid as long as the neutron-separation energy is larger than Qp. 
Since, at excitation energies higher than Sn, de-excitation through neutron emission is favored over decay by y-ray 
emission, thus implying Ftot — Fn ^ T-yjFf [371 148 ■ For all nuclei mentioned in Table however, Qp is below Sn- 

13), can thus also be derived for neutron-rich /3DF precursors by 


An approximation of 7i/2p,/3DF) similar to equation 
taking into account above considerations. However, considering the limited experimental information on /3DF in the 
neutron-rich region, a detailed derivation is omitted in this paper. 
















Figure 4. (Color online) Plots of ri/ 2 p,/ 3 DF versus [Qp — B{) for different fission models as listed in Table|^ The closed symbols, 
representing reliable values for Ti/ 2 p,^DF in Table 11 are used for a linear fit with equal weights to the data points. Other data 
from Table [I] are indicated by the open symbols. The color code represents the different regions of the nuclear chart for which 
/?DF has been experimentally observed : the neutron-deficient lead region (red), the neutron-deficient (black) and neutron-rich 
(blue) uranium region. 


IV. CONCLUSIONS 


Recent experiments have measured the /3DF of 9 precursor nuclei in the neutron-deficient lead region. Because 
of the dominant a-decay branch in most of these nuclei, /3DF probabilities are extracted with large experimental 
uncertainties. In contrast, the partial half-life for /3DF can be determined with a better accuracy. In addition, 
’^i/ 2 p,/ 3 DF can be easily derived from the /3DF probability by using equation (0. 

A systematical evaluation of /3DF partial half-lives was performed by using hssion barriers deduced from four different 
models for a broad range of nuclei in the lead and uranium regions. A linear relation between logic(T i/ 2 p,/ 3 df) ^^nd 
iQ /3 ~ Bf) was observed for neutron-deficient precursor nuclei, when using tabulated fission barriers from the TF 
or LDM approach, and to a lesser extent for FRLDM and ETFSI barriers. This linear trend persists for valnes of 
2 ^i/2p,/3DF spanning over 7 orders of magnitude and a wide variety of precursor nuclei going from i'^®Tl to ^^®Es with 
N/Z ratios of 1.20 and 1.51 respectively. This observation may help to assess /3DF branching-ratios in very neutron- 
rich nuclei, which are inaccessible by present experimental techniques but might play a role in the fission-recycling 
mechanism of the r-process nucleosynthesis. 
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